WILLIAMS' DECOMPOSITION OF THE LEVY CONTINUUM RANDOM 
TREE AND SIMULTANEOUS EXTINCTION PROBABILITY FOR 
POPULATIONS WITH NEUTRAL MUTATIONS 



ROMAIN ABRAHAM AND JEAN-FRANgOIS DELMAS 

Abstract. We consider an initial Eve-population and a population of neutral mutants, 
such that the total population dies out in finite time. We describe the evolution of the 
Eve-population and the total population with continuous state branching processes, and the 
neutral mutation procedure can be seen as an immigration process with intensity propor- 
tional to the size of the population. First we establish a Williams' decomposition of the 
genealogy of the total population given by a continuum random tree, according to the an- 
cestral lineage of the last individual alive. This allows us to give a closed formula for the 
probability of simultaneous extinction of the Eve-population and the total population. 



1. Introduction 

We consider an initial Eve-population whose size evolves as a continuous state branching 
process (CB), = {Y^^,t > 0), with branching mechanism V'Eve- We assume this population 
gives birth to a population of irreversible mutants. The new mutants population can be 
seen as an immigration process with rate proportional to the size of the Eve-population. We 
assume the mutations are neutral, so that this second population evolves according to the 
same branching mechanism as the Eve-population. This population of mutants gives birth 
also to a population of other irreversible mutants, with rate proportional to its size, and so 
on. In [2], we proved that the distribution of the total population size Y = (Yf,t > 0), which 
is a CB with immigration (CBI) proportional to its own size, is in fact a CB, whose branching 
mechanism ^p depends on the immigration intensity. The joint law of {Y^, Y) is characterized 
by its Laplace transform, see Section 14.1.41 This model can also be viewed as a special case 
of multitype CB, with two types and 1, the individuals of type giving birth to offsprings 
of type or 1, whereas individuals of type 1 only have type 1 offsprings, see [13\ |6] for recent 
related works. 

In the particular case of Y being a sub-critical or critical CB with quadratic branching 
mechanism {^p{u) = oqu + (3v?^ /9 > 0, ao ^ 0)> the probability for the Eve-population to 
disappear at the same time as the whole population is known, see [T7] for the critical case, 
ao = 0, or Section 5 in [2] for the sub-critical case, ao > 0. Our aim is to extend those results 
for the large class of CB with unbounded total variation and a.s. extinction. Formulas given 
in [2] could certainly be extended to a general branching mechanism, but first computations 
seem to be rather involved. 

In fact, to compute those quantities, we choose here to rely on the description of the 
genealogy of sub-critical or critical CB introduced by Le Gall and Le Jan [12j and developed 
later by Duquesne and Le Gall [7], see also Lambert [TO] for the genealogy of CBI with 
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constant immigration rate. Le Gall and Le Jan defined via a Levy process X the so-called 
height process H = {Ht,t > 0) which codes a continuum random tree (CRT) that describes 
the genealogy of the CB (see the next section for the definition of H and the coding of the 
CRT). Initially, the CRT was introduced by Aldous [4] in the quadratic case: ip{X) = A^. 
Except in this quadratic case, the height process H is not Markov and so is difficult to 
handle. That is why they also introduce a measure- valued Markov process {pt,t > 0) called 
the exploration process and such that the closed support of the measure pt is [0, Ht] (see also 
the next section for the definition of the exploration process). 

We shall be interested in the case where a.s. the extinction of the whole population holds 
in finite time. The branching mechanism of the total population, Y, is given by: for A > 0, 

(1) V(A) = aoA + /3A2 + / Tr(di) (e-^^ -1 + Xi) , 

J{0,oo) ^ ' 

where ao ^ 0, /3 > and vr is a Radon measure on (0,oo) such that /(q oo)(^ ^('^^) 

< oo. 

We shall assume that Y is of infinite variation, that is /? > or i'iT{di) = oo. We shall 
assume that a.s. the extinction of Y in finite time holds, that is, see Corollary 1.4.2 in [T], 
we assume that 

(2) / — - < oo. 



We suppose that the process Y is the canonical process on the Skorokhod space B(M+, M+) of 
cadlag paths and that the pair (Y, Y^) is the canonical process on the space D(M_|_, M_|_)^. Let 
denote the law of the pair (y,y°) (see [2]) started at (10,^0°) = (^'^)- The probability 
measure P^; is infinitely divisible and hence admits a canonical measure N: it is a cj-finite 
measure on B(M_|_,M+)^ such that 

(y,yO)(J^(r,yO'*) 

where ((y*, y*^'*), i G /) are the atoms of a Poisson measure on D(M_|_, M_|_)^ with intensity 
xN^dY, dY^). In particular, we have 

(3) E^[e-^^'] = exp(-rEN[l - e"^^*]) 
and u{X,t) = N[l - e"-^"*^*] is the unique non-negative solution of 

(4) / -— = t, for t > and A > 0. 

Juix,t) i^{v) 

Let Ty = inf {t > 0; Yj = 0} be the extinction time of Y. Letting A go to co in the previous 
equalities leads to 

Px{ty <t) = exp -xN[Ty > t], 
where the positive function c{t) = N[ry > t] solves 

r°° dv 

(5) / 77^=*, fort>0. 

Jc{t) ■^[V] 

Let us consider the exploration process {pt,t > 0) associated with this CB. We denote by 
N its excursion measure. Recall that the closed support of the measure pt is [0,i7f], where 
H is the height process. Let be the total local time at level a of the height process H 
(well-defined under N). Then, the process {L"',a > 0) under N has the same distribution as 
the CB Y under N. 



WILLIAMS' DECOMPOSITION 



3 



We decompose the exploration process, under the excursion measure, according to the 
maximum of the height process. In terms of the CRT, this means that we consider the 
longest rooted branch of the CRT and describe how the different subtrees are grafted along 
that branch, see Theorem 13.31 When the branching mechanism is quadratic, the height 
process H is a Brownian excursion and the exploration process pt is, up to a constant, the 
Lebesgue measure on In that case, this decomposition corresponds to Williams' 

original decomposition of the Brownian excursion (see llSj). This kind of tree decomposition 
with respect to a particular branch (or a particular subtree) is not new, let us cite [HI [Uj for 
instance, or [16l [151 [H] for related works on superprocesses. 

We present in the introduction a Poisson decomposition for the CB only, and we refer to 
Theorem l3.3l for the decomposition of the exploration process. Conditionally on the extinction 
time Ty equal to m, we can represent the process Y as the sum of the descendants of the 
ancestors of the last individual alive. More precisely, let M'{dl,dt) = 6^£.^t.){d£, dt) be a 

iei 

Poisson point measure with intensity 

lM(i)e-'^(™-*) £7T{di)dt, 

and 

(6) (dt) = Y,^i^tMt) + 2pl[o^m)it)dt. 

Let Nt{dY) denote the law of {Y(s — t),s > t) under N and ^^'^(tj,yj) be, conditionally 
on J\f', a Poisson point measure with intensity 

{dt)Nt[dY,l{^<^y] 

where Nt[dY, 1|^<„}] denotes the restriction of the measure Nf to the event {ry < m}. 
The next result is a direct consequence of Theorem 13.31 

Proposition 1.1. The process Yljej^'^ distributed as Y under 1^, conditionally on {ry = 
m}. 

Let Tyo = inf{t > 0;!^" = 0} be the extinction time of the Eve-population. In the 
particular case where the branching mechanism of the Eve-population is given by a shift of 
^|;■. 

(7) i;E.e{-)=^PiO + -)-m, 

for some 6 > and /3 = 0, the pruning procedure developed in [l] gives that the nodes of 
width ii correspond to a mutation with probability 1 — e~^^\ As /? = there is no mutation 
on the skeleton of the CRT outside the nodes. In particular, we see simultaneous extinction 
of the whole population and the Eve-population if there is no mutation on the nodes in 
the ancestral lineage of the last individual alive. This happens, conditionally on Kmax, with 
probability 
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Integrating w.r.t. the law of J\f gives that the probabihty of simultaneous extinction, condi- 
tionally on {ty = m}, is under N, given by 

N[ryo = m\TY = m] = exp - y l[o,„)(t) e'^'^'"^*) l'K{dl)dt 1 - e^' 

= exp - / [V''(c(m -t)+e)- 7p'{c{m - t))] dt 
Jo 

= exp — / (f){c{t)) dt, 



where (p = V'Eve — V'- Now, using that the distribution of {Y^,Y) is infinitely divisible 
with canonical measure N, standard computations for Poisson measure yield that PxiTyo = 
m\TY = m) = N[ryo = m|ry = m] that is 

I'm 

Px{ty" = w-lry = m) = exp— / 4''{c{t)) dt. 

Jo 

Notice that this formula is also valid for the quadratic branching mechanism {ip{u) = oqu + 
/3m^ /3 > 0, ao > 0), see Remark 5.3 in [2]. 

In fact this formula is true in a general framework. Following [2] , we consider the branching 
mechanisms of the total population and Eve-population are given by 

2 , f _^j/)\r„-A£ 



ipiX) = aoX + PX^+ / 7T{de)[e-^' -1 + X£], 



'0Eve(A) = OEveX + fSX^ + / 7rEvc(d^) [e"^^ "1 + A^], 

J(0,oo) 

and the immigration function 



<PiX) = VEve(A) - V'(A) = ai^mA + / - e ), 

J{0,od) 

where aimm = OEve — oo — /(g oo) ^^(^^) ^ ^ = ^Eve + i^, where ttevg and v are 
Radon measures on (0, oo) with J^^ lu{dl) < oo. Notice the condition J^q lu{dtj < oo is 
stronger than the usual condition on the immigration measure, J^g oo)^''' ^ ^) i^{dl) < oo, but 
is implied by the requirement that J^-^ £v{dtj < J^^ i7r{d£) < oo. 

Inspired by Theorem 13.31 we consider M{dl,dt,dz) = ''^^6(^£^^t^^Zi){d£,dt,dz) a Poisson 
point measure with intensity 

(8) i[o,m){t) e-'^""-'^ I [7rEve(d^)5o(d^) + v{dl)bx(dz)\ dt. 

Intuitively, the mark Zi indicates if the ancestor (of the last individual alive) alive at time ti 
had a new mutation (zj = 1) or not (zi = 0). Note however that if /3 > we have to take 
into account mutation on the skeleton. More precisely, let Ti = min{tj,2:j = 1} be the first 
mutation on the nodes in the ancestral lineage of the last individual alive and let T2 be an 
exponential random time with parameter aimm independent of M. The time T2 corresponds 
to the first mutation on the skeleton for the ancestral lineage of the last individual alive. We 
set 

^g-^ jTo = min(ri,r2) if min(Ti, r2) < m, 

1 To = -|-oo otherwise. 

In particular there is simultaneous extinction if and only if Tq = +00. 
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For t > 0, let us denote by Nt{dY^, dY) the joint law of {{Y^{s - t),Y{s - 1)), s>t) under 
N. Recall Kmax given by Conditionally on M and T2, let '^{tj,YO,j,-Kj) be a Poisson 

point measure, with intensity 
We set 

(10) {Y'\y') = J2 {Y''',Y^) + E (0'^')- 

tj<To tj>To 

We write Qm for the law of (Y'^,Y') computed for a given value of m. 

Theorem 1.2. Under Qm, {Y'^,Y') is distributed as {Y^,Y) under N[-\ty = m] , or equiv- 
alently, under j^"^ \c' {m)\Qm{')d'm, iY' ,Y') is distributed as {Y^,Y) under N. 

Let us remark that this Theorem is very close to Theorem 13.31 but only deals with CB and 
does not specify the underlying genealogical structure. This is the purpose of a forthcoming 
paper [3j where the genealogy of multi-type CB is described. 

Intuitively, conditionally on the last individual alive being at time m, until the first mu- 
tation in the ancestral lineage (that is for tj < Tq) , its ancestors give birth to a population 
with initial Eve type which has to die before time m, and after the first mutation on the 
ancestral lineage (that is for tj > Tq), there is no Eve-population in the descendants which 
still have to die before time m. 

Now, using that the distribution of {Y^,Y) is infinitely divisible with canonical measure 
N, standard computations for Poisson measure yield that Px{tyo = m\TY = ra) = N[ryo = 
m|ry = m]. As 

N[TyO = m\TY = m]= Qm{To = +00) 

= Qm{Ti = +oo)Q^(r2 > m) 

we deduce the following Corollary. 

Corollary 1.3 (Probability of simultaneous extinction). We have for almost every m > 

Pa;(Tyo = m\TY = m) = exp — / (j)' {c{t)) dt, 

Jo 

where c is the unique (non-negative) solution of (0). 

The paper is organized as follows. In Section O we recall some facts on the genealogy 
of the CRT associated with a Levy process. We prove a Williams' decomposition for the 
exploration process associated with the CRT in Section [3l We prove Theorem 11.21 in Section 
m Notice that Proposition 11.11 is a direct consequence of Theorem 11.21 

2. Notations 

We recall here the construction of the Levy continuum random tree (CRT) introduced 
in [12\ ITT] and developed later in [7]. We will emphasize on the height process and the 
exploration process which are the key tools to handle this tree. The results of this section 
are mainly extracted from [7]. 
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2.1. The underlying Levy process. We consider a R-valued Levy process {Xt,t > 0) with 
Laplace exponent ^ (for A > E [e"^-^*] = e*'''^^)) satisfying ([I]) and ([2]). Let / = {It,t > 0) 
be the infimum process of X, It = mio<s<t Xg, and let S = {St,t > 0) be the supremum 
process, St = supo<s<t^s- We will also consider for every < s < t the infimum of X over 
[s,t]: 

If = inf Xr. 

s<r<t 

The point is regular for the Markov process X — I, and — / is the local time of X — / at 
(see [5], chap. VII). Let N be the associated excursion measure of the process X — I away 
from 0, and a = inf{t > 0; Xt — It = 0} the length of the excursion of X — I under N. We 
will assume that under N, Xq = Iq = 0. 

Since X is of infinite variation, is also regular for the Markov process S — X. The local 
time, L = {Lt,t > 0), of S" — X at will be normalized so that 

where = inf{s > 0; > t} (see also Theorem VII. 4 (ii)). 

2.2. The height process and the Levy CRT. For each t > 0, we consider the reversed 
process at time t, = (xi*\ < s < t) by: 

Xi'^ =Xt-X^t-s)- if 0<s<t, 

and X^^*^ = Xt- The two processes {Xs^\o < s < t) and (^^,0 < s < t) have the same law. 
Let be the supremum process of and L^*) be the local time at of 5^*) - X^*) with 
the same normalization as L. 

Definition 2.1 (^, Definition 1.2.1 and Theorem 1.4.3). There exists a process H = {IIt,t > 
0), called the height process, such that for all t > 0, a.s. Ht = Lf\ and Hq = 0. Because of 
hypothesis the height process H is continuous. 

The height process {IIt,t € [0, o"]) under N codes a continuous genealogical structure, the 
Levy CRT, via the following procedure. 

(i) To each t G [0, a] corresponds a vertex at generation Ht- 

(ii) Vertex t is an ancestor of vertex t' if Ht = H^t,t'] ) where 

(11) Hyt^t'] = inf{^„, ue[tM\t\J t']]. 

In general II[t,t'] is the generation of the last common ancestor to t and t' . 

(iii) We put d{t, t'j = Ht + Hf - 2H)^t,t'] and identify t and t' {t ~ t') if d{t, t') = 0. 

The Levy CRT coded by H is then the quotient set [0,cr]/ equipped with the distance d 
and the genealogical relation specified in (ii). 

Let (ts, s > 0) be the right continuous inverse of — /: = inf{t > 0; —It > s}. Recall that 
— / is the local time of X — / at 0. Let denote the local time at level a oi H until time i, 
see Section 1.3 in [^Tj. 

Theorem 2.2 (^, Theorem 1.4.1). The process (L'^^,a > 0) is under ¥ (resp. defined 
as Y under Fx (resp. NJ. 

In what follows, we will use the notation N instead of N for the excursion measure to stress 
that we consider the genealogical structure of the branching process. 
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2.3. The exploration process. The height process is not Markov. But it is a simple 
function of a measure-valued Markov process, the so-called exploration process. 

If E' is a Polish space, let B{E) (resp. B^{E)) be the set of real-valued measurable (resp. 
and non-negative) functions defined on E endowed with its Borel fi-field, and let Ai{E) (resp. 
A4f{E)) be the set of cr-finite (resp. finite) measures on E, endowed with the topology of 
vague (resp. weak) convergence. For any measure fi S M{E) and / € B^{E), we write 

(/.,/) = J f{x)f,{dx). 

The exploration process p = {pt,t > 0) is a 7V4j(M+)-valued process defined as follows: for 
every / G S+(M+), 

{ptJ)= [ dsl!f{Hs), 

J[o,t] 

or equivalent ly 

(12) pt{dr)= ^ m - Xs-)5HMr) + f31[o,H,]{r)dr. 

0<s<t 

In particular, the total mass of pt is {pt, 1) = Xt — It- 
For p G A1(]R+), we set 

(13) H{p) = sup Supp p, 

where Supp p is the closed support of p, with the convention H{0) = 0. We have 

Proposition 2.3 ([7j, Lemma 1.2.2). Almost surely, for every t > 0, 
. H{pt) = Ht, 

• Pt = if and only if Ht = 0, 

• if pt^ 0, then Supp pt = [0,Ht]. 

In the definition of the exploration process, as X starts from 0, we have po = a.s. To state 
the Markov property of p, we must first define the process p started at any initial measure 
p G Mf(R+). 

For a £ [0, {p, 1)], we define the erased measure kaP by 

kap{[0, r]) = p{[0, r]) A {{p, 1) - a), for r > 0. 

If a > {p, 1), we set kaP = 0. In other words, the measure kaP is the measure p erased by a 
mass a from the top of [0, H{p)]. 

For v,p £ A^j(M+), and p with compact support, we define the concatenation [piv] G 
M.f(R+) of the two measures by: 

{[pMJ) = (^,/> + {yJ{H{p) + •)>, / e 

Finally, we set for every p G A1j(]R+) and every t > = \k^i^p, pt]. We say that 

(Pt'ji > 0) is the process p started at p^ = p, and write for its law. Unless there is an 
ambiguity, we shall write pt for p^. 

Proposition 2.4 (|7], Proposition 1.2.3). The process {pt,t >0) is a cad-lag strong Markov 
process in 7V4j(M+). 

Notice that N is also the excursion measure of the process p away from 0, and that a, the 
length of the excursion, is N-a.e. equal to inf{t > 0; pt = 0} . 
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2.4. The dual process and representation formula. We shall need the A^j(M+)- valued 
process r] = {r]t,t > 0) defined by 

rhidr)= iXs-inSHMr)+PMoM^)dr. 

0<a<t 
Xs-<It 

The process r/ is the dual process of p under N thanks to the following time reversal property: 
recall a denotes the length of the excursion under N. 

Proposition 2.5 ([7j, Corollary 3.1.6). The processes {{ps,r]s);s > 0) and ((f?(o— s)-; P{(t-s)-); 

s > 0) have the same distribution under N. 

It also enjoys the snake property: for all t > 0, s > 

(/5t,^t)[0,H[,,,]) = {ps,Vs)[0,Hi,^,]), 

that is the measures p and r] between two instants coincide up to the minimum of the height 
process between those two instants. 

We recall the Poisson representation of {p, rj) under N. Let Af^{dx d£ du) be a Poisson point 
measure on [0, +cx3)^ with intensity 

dx £7r{d£)l[Q^ij {u)du. 

For every a > 0, let us denote by Mq the law of the pair {pa-,Va) of finite measures on 
defined by: for / G e+(R+) 



{PaJ) = J J\f*{dxd£du)l[Q^a]ix)u£f{x), 
{uaj) = [ K{dxd£du)^o^,^{x)£{l-u)f{x). 



We finally set M = /+°° da e-'^^" M^. 

Proposition 2.6 ([.7j, Proposition 3.1.3). For every non-negative measurable function F on 



N 



M{dpdiy)F{p,u) 



F{pt,rjt)dt -- 

where a = inf{s > 0; ps = 0} denotes the length of the excursion. 

We can then deduce the following Proposition. 
Proposition 2.7. For every non-negative measurable function F on A4f{ 

g-«oa f Maidpdi^)F{p,u), 



N / Fipt,m)dL'} 

Jo J J 

where a = mf{s > 0; ps = 0} denotes the length of the excursion. 

3. Williams' decomposition 

We work under the excursion measure. As the height process is continuous, its supremum 
-H'max = sup{Hr;r € [0,cr]} is attained. Let Tmax = inf{s >0;Hs = -H'max}- 

For every m > 0, we set Tm{p) = inf{s > 0, Hs{p) = m} the first hitting time of m for the 
height process. When there is no need to stress the dependence in p, we shall write Tm for 
Tm (p) • Recall the function c defined by ([5]) is equal to 

(14) c(m) = N[r^ < 00] = N[i7„,ax < m]. 
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We set pd = {pT^,^+s, s > 0) and Pg = {P(t^,^-s)+, « > 0), where x+ = max(x, 0). 
For every finite measure with compact support p, we write P* for the law of the exploration 
process p starting at p and killed when it first reaches 0. We also set 

P;^ := hm P;^( . I H{p) < H^,,, < Hip) + e). 

We now describe the probability measure P* via a Poisson decomposition. Let (Qj,/5j), i € / 
be the excursion intervals of the process X — I away from (well defined under P* or under 
P*). For every i £ I, we define hi = Ha. and the measure- valued process by the formula 



{ptJ) = f{x- hi)p(^a,+t)ApMx)- 

J (hi,+oo) 

We then have the following result. 

Lemma 3.1. Under the probability P*, the point measure ^(hi,p'-) Poisson point 

iei 

measure with intensity p{dr)N[-, H^ax <m — r]. 

Proof. We know (cf Lemma 4.2.4 of [7j) that the point measure ^^'^(hj.p*) is under P* a 

Poisson point measure with intensity p{dr)N{dp). The result follows then easily from standard 
results on Poisson point measures. □ 

Remark 3.2. Lemma |3 . 1 1 gives also that, for every finite measure with compact support p, if 
we write pa = p{- H [0, a]), 

^1= hm f;j,\h^,,,<h{p)). 

Theorem 3.3 (Williams' Decomposition). 

(i) The law of H^ax is characterized by N[Hmax < rn] = c{m), where c is the unique 
non-negative solution of 

(ii) Conditionally on H^ax = "i, the law of {pTmaxi'HTmax) is under N the law of 

I "^Viridt^ +/31[o,m](t)dt,^(l -v.i)ri5t^ + /?l[o,m] (i)di | , 
Kiel iei J 

where ^^(^^.^r.^ti) is a Poisson measure with intensity 

l[o,i]{v)lioM(.^)e-''^"'"'Uvr7r{dr)dt. 

(iii) Under 'N, conditionally on H^ax = fn, and {pTmax^'HTmax)^ {Pd, Pg) o.re independent and 
Pd fresp. Pa) is distributed as p (resp. n) under ¥* fresp. P* ). 

y rimax Umax 

Notice (i) is a consequence of (|14p . Point (ii) is reminiscent of Theorem 4.6.2 in [7] which 
gives the description of the exploration process at a first hitting time of the Levy snake. 

The end of this section is devoted to the proof of (ii) and (iii) of this Theorem. 

Let m > a > be fixed. Let e > 0. Recall Tm = mf{t > 0;Ht = m} is the first hitting 
time of m for the height process, and set Lm = sup{i < a;Ht = m} for the last hitting time 
of m, with the convention that inf0 = +00 and sup0 = +00. We consider the minimum of 
H between Tm and L^: -f^[T™,L^] = ram{Ht]t € [Tm,Lm]}. 
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We set = (/>T^^,,+t,t > 0), with 
the path of the exploration process on the right of Tmax after the hitting time of a, and 



P 



(9) 



(/'(Lmax.a-t)-'* ^ 0), with Lmax.a = sup{t < T^aKx] Ht = a}, the returned path of the 
exploration process on the left of Tmax before its last hitting time of a. Let us note that, by 
time reversal (see Proposition I2.5p . the process p^^^ is of the same type as rj. This remark 
will be used later. 

To prove the Theorem, we shall compute 



N 



,a])-'-{m<_fi'max<m,+e} 



and let e go down to 0. We shall see in Lemma [3^ that adding Iro- ^ in the integrand 
does not change the asymptotic behavior as e goes down to 0. Intuitively, if the maximum of 
the height process is between m and m + e, outside a set of small measure, the height process 
does not reach level a between the first and last hitting time of m. So that we shall compute 
first 



(15) 



A = N 



Fl (p(^) ) F2 ip^^^)F3 [pT^^^ I [0 )F4 (7?T„,ax I [0,a] )'^{H[t^ >a,m<H„,ax<»n+e} 



Notice that on {HyTm,Lm\ > ^^'^^ ^max,a = Tm,a ■= inf{s > Tm,Hs{p) = a} and, 

from the snake property, />T,„ax | [o,a] = PT™|[o,a] and ??Tn,ax|[o,a] = ^T™|[o,a]' so that 



A = N 



Fi{/a^)F2i{pT^^^+t,t> 0))F3(/5T^|[o,,])F4(r?I 



Let US remark that, we have 



-'-{^^[Tm,Lml>«i™<^^max<"!.+e} ^ {m<snp{Hu ,0<u<Tm,a}<m+e}^ {snp{Hn ,U>Tm,a} <m} ■ 

By using the strong Markov property of the exploration process at time Tm,a; we get 



P'^iP^^^)P4.ir]Trn\[0,a])'^{m<sup{H^,0<u<Trn,a}<m+s}F3{PTrn\[0,^^^^ .[^2(p)l{//^ax<m}] 



and so, by conditioning, we get 



^ = N 



Fl (p(^) ) Fi (7?T„ I [0^^] ) G2 {PT„, I [o,a] ) 1{/^[T„ , 



l^^ >a,m<Hmax<m+e} 



where G2{p) = -F3(p)E* [F2(p)|ffmax < ""t-]- Using time reversibility (see Proposition 12.5 
and the strong Markov property at time T^.a again, we have 



A = n 

= N 



i^l(p('))F4(pT„|[o,a])G'2(??T„ 

I [0,a] ' {H[Trn,Lm,] >a,m<Hmayc<m+e} 
Gl [PTm I [0,a] ) G2 {riTrr, | [0,a] ) 1{J^[t™ ,L„] >a,m<//max<m+e} 



where Gi{p) = Fi{p)W^[Fi{p)\H^^^ < m]. 

Now, we use ideas from the proof of Theorem 4.6.2 of [7|. Let us recall the excursion decom- 



position of the exploration process above level a. We set r° = inf < r, / dul 



> S 



Let £a be the fi-field generated by the process {ps, s > 0) := {pr^, s > 0). We also set 
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Let {ai, Pi), i e I he the excursion intervals of H above level a. For every i e I we define 
the measure-valued process by setting 

f (Ps> = /(a,+oo) Pai+sidr)(p{r -a) if < s < A - "i, 

1 = if s = or s > /3i — aj, 

and the process 77* similarly. We also define the local time at the beginning of excursion 
by £i = Lq.. Then, under N, conditionally on £a, the point measure 

is a Poisson measure with intensity lyQ^La^{t)d(M\dp drj\. 
In particular, we have 



A = n 



En 

^{T^ipi )=+oo} 



GliPai)G2i'nai)'^{m<H, 



i(p»)<m+e} 



Let us denote by (r^,^ > 0) the right-continuous inverse of {Lg,s > 0). Palm formula for 
Poisson point measures yields 



A = N 



= N 



N 



{m<i?max (p*)<m+£} 



diGi{prf)G2{rirf)Hm < H^^ <m + e]n 



A time-change then gives 



(16) 



A = v{m — a,s)N 



Jo 



where v{x, e) = c{x) — c{x + e) = N[x < -ffmax < x + e]. We have 



A = v{m — a,£)N 
= v{m — a,£)N 



r dL1G^{ps)G2{r]s)^ 
Jo 



-c{m-a)L'i ^-c{m-a){L%-L'i) 



dL^Gi(p,)G2(%)e 



-c{m-a)L'^^ -(ps,N[l-e-'=("»-°)-^S ']) 



where we used for the last equality that the predictable projection of e "^^^^ is given by 
g-(ps,N[i-e ^^CT ]) jsjQ^jgg \jy using the excursion decomposition above level < r < m, 
we have 

c(m) = N[r„ < 00] = N[l - 

In particular, we get 



A = v{m — a,£)N 
Using time reversibility, we have 
A = v{m — a,£)N 



dL«Gi(p,)G2(r?,)e 



r dLlGi{'ns)G2{ps 
Jo 



-c(m-a)(L--L-) (r,,,c(m-)> 
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Similar computations as those previously done give 

dL«Gi {ris ) G2 {ps ) e- 







dL:Gi{ps)G2ivs)e 



-{ps+Vs,c{m--)) 



A = v{m — a, e)N 

= v{m — a, e)N 

Using Proposition 12.7^ we get 

A = v{m- a, e) e"""'^ J Ma{dpdiy)Gi{n)G2{u) e-</^+^''=("^-)) . 

We can give a first consequence of the previous computation. 
Lemma 3.4. We have 

c(m — a) — dm — a + e) 



N[i7[T^_i^] > a,m< Hmax <m + e]= c{m)- 



d{m — a) 



Proof. Taking Fi = F2 = F3 = F4 = 1 in ((H]), we deduce that 
N[iT[j;^_i^] > a, m < i^max <m + e] = v{m - a, e)N 



j^a ^-c{m-a)L'^ 



Let ao > and let us compute B{ao,a) 
that 



. Thanks to Theorem 12.21 notice 



B{ao,a) = N 
On the other hand, we have 



9aoN[l-e 



-c{ao)Ya 



c'(ao) 



c(a + ao) = Wa+ao > 0] = N[l - Ey„[y,„ = 0]] = N 



1 - e 



-Yac{ao) 



where we used the Markov property of Y at time a under N for the second equality and ( 

with A going to infinity for the last. Thus, we get B(aQ,a) = — ^-77— r-^- We deduce that 

c'(ao) 

HH[Tm,L,„] > a,m< i?max < m + e] = v{m - a, e)B{a - m, a) 

. dm — a) — dm — a + e) 

= c (m) ^ . 

c'(m — a) 



□ 



Since ^1,^2,^3 and F4 are bounded, say by C, we have |A — ^ C^^[H[Tm,Lm] < CL,m < 
Hma.x < m + e]. From Lemma 13. 4[ we deduce that 



r l^-^ol 

iim 

e^O N[m < -ffmax < m + e] 



1 — lim 



N[-H'[T^,L^] > Q, m < ffmax < m + e\ 
N[m < i^max <m + e\ 
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We deduce that 



N 

lim 



I [0,a] J -'-{rri<_ffmax<m+e} 



c'(m) 7 
_ / Ma{dfi du)Gi {fi)G2{v) e-{^^+'^^<rn-)) 



J Maidfi du) e-(/'+''.c{'»-)) 

la{d^ldv)Gl{^i)G2{v) 

K{dpdu)Fi{ij)El[F,{p^''^)\H„,,, < m]F3(//)E;;[F2(p('^))|if„,ax < m], 



where 



p{dt) = "^Uilidt^ + pi[Q^ci]i't)dt 

v{dt) = Y,{'^ - Ui)^i5u+l31yQ^a]{t)dt, 
iei 

and Yliei ^{xi,et,ti) is under Mq a Poisson point measure on [0, +00)^ with intensity 

l[o,a]{t)dtie-''^"'-'^7T{di)l[o,,^iu)du. 

Standard results on measure decomposition imply there exists a regular version of the 
probability measure N[ • |i^max = w-] and that, for almost every non-negative m, 

N[ • |i?max = m] = lim N[ • |m < F^ax <m + e]. 

This gives (ii) and (iii) of Theorem l3.3l since Fi,F2, F-^^F^ are arbitrary continuous functionals 
and by Remark [ 



4. Proof of Theorem 11.21 

The proof of this Theorem relies on the computation of the Laplace transform for {Y'^ ^ Y') 
and is given in the next three paragraphs. The next paragraph gives some preliminary 
computations. 

4.1. Preliminary computations. 

4.1.1. Law oJTq. Recall the definition of Qm as the law of (y'°, ¥') defined by ^ and Tq 
defined by Q as the first mutation undergone by the last individual alive. 
For r < m, we have 

Qmin e[r,r + dr],To = T2) = Qm{T2 G [r, r + dr])Q„,{Ti > r) 

= dr ai^^ e-"i--^ exp - [ dt [ q-^^^^*) eu{d£) 

Jo J{0,oo) 
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and, with the notation (po{X) = (j){X) — aimm'^! 
Qm{To e[r,r + dr],To = Ti) = Q^(r2 > r)Q„(ri G [r, r + dr]) 

= dr (p'oicim - r)) e-"!--"^ exp - T dt [ ^-H^-t) i,y(^di) 

Jo J(0,oo) 

= dr (p'o{c{m - r)) e" /o -^'W"^-*)) '^^ . 
In particular, we have for r < m 

Qm{To G [r,r + dr]) = dr (/.'(c(m - r)) e" i'o '^'W'""*)) '^K 

and 

(17) Q,„(ro > r) = e" -fo <^'(c(™-t)) dt _ 



Notice we have Qm(2o = oo) = exp— / (j)'{c{t)) dt. 

Jo 



4.1.2. Conditional law of M given Tq. Recall M is under a Poisson point measure with 
intensity given by ([8|). Conditionally on {Tq = r, Tq = T2}, with m > r > 0, TV is under Qm 
a point Poisson measure with intensity 

l[o,r)(i) e-^'^^'"-*) ^^^,e{d^)8o{dz)dt+ 

l{r,m)(i) e-^'=('"-*) ^ [7rEve(d^)5o(d^) + v{dt)5i{dz)] dt. 

Conditionally on {Tq = r, Tq = Ti}, with r < m, TV is distributed under Qm as TV + (5(j^^^,^i) 
where TV is a point Poisson measure with intensity 

l[o,r)(t) e-^'^^'"-*) £7rEve((i^)<5o(d^)di 

+ l(r,-m){t) e-^'=("'-*) ^ [7rEvc(d^)<Jo((i-z) + v{dt)5i{dz)] dt, 
and L is a random variable independent of TV with distribution 

Conditionally on {Tq = cxo}, TV is under Qm a point Poisson measure with intensity 

l[o,m)(i) e~^^(™'*) ^7rEve(d^)5o(d^)di. 

4.1.3. Formulas. The following two formulas are straightforward: for all x,7 > 0, 

(18) VEve(^ + 7) - V-EvelT) = 2(3x + [ e"^^ ^7rEve(d^)[l " 6"^^, 

J(0,oo) 

(19) V'(a; + 7) -^'(7) = 2/32;+ / e-^^£7r(d^)[l-e-'^^], 

7(0,00) 

Finally we deduce from ([5]) that ■i/'(c) = — c', tp'{c)c' = —c" and 

(20) / V''(c) = -log(c'). 



WILLIAMS' DECOMPOSITION 



15 



4.1.4. Laplace transform. Recall Ty = ini{t > 0; 1^ = 0} is the extinction time of Y. Let 
jUEve and //Total be two finite measures with support a subset of a finite set A = {ai, . . . , a„} 
with = ao < fli < • • • < On < ^n+i = oo. For m G (0, +oo) \ A, we consider 

= N[l - e- Z^-' '^-.a.W 

By noticing that N-a.e. l{ry<m-t} = l™A^oo exp — / y^._^ fi^{dr)^ where fi^{dr) = \5m{dr)^ 
we deduce from Lemma 3.1 in [2j that iwYa") 

w^) are right continuous and are the unique 
non-negative solutions of : for G {0, . . . , n}, m € (a^, a^+i), t E (—00, m), 



(21) w*^{t)+ I ^{w*^{r))dr = I fJ-Totid{dr) + c{m - ak), 



(22) Wm{t) + / ■ipEve{wm{r))dr 

J[t,ak] 

= / //Eve('^?^) + / /UTotal ((ir) + c(m - afc) + / (j){w^{r))dr. 

J[t,a^] J[t,ak\ -^[t.^fc] 

We define 

(23) dm = maxjafc; < m, e {0, . . . , n}}. 
Notice that Wm{t) = w^it) = c{m — t) for t S {dm,m). 

4.2. Proof of Theorem [ll2l 

4.2.1. 74im. Theorem 11.21 will be proved as soon as we check that the following equality 
w{0) = / -c'(m)Q„[l - e-^^'- /^EveK)-/ y; /.xotaiW]^^ 



JO 

holds for all the possible choices of measures /^Eve and //Total satisfying the assumptions of 
Section 14.1.41 with w = w^o defined by (p2|) . 

Notice the integrand of the right-hand side is null for m < ai . Let A denote the right-hand 
side. We have for < e < oi: 

/■oo 

A = / dm (-c'(m))Q„[l - e"/^'-- /^evc (*•)-/ y; MTotai(<ir)j 

/>oo 

= c(e)+ / dmlA-{m)c'{m)Qm[Z], 



with, thanks to the definition (161) of k,, 



Z = exp - / Kmax(di) [ntl{f<To} + nll{t>To}] 

Jo 



and 



N[(l - e-/^-' '^-°-'("'-))l{^<„.-o] = - - t), 



with {wm,Wm) the non- negative solutions of ([2T]) and (p2|) . Notice that Wm{t) = w^it) 
c{m — t) for t G (a^, w-) and thus nj = = when t G {dm, m). 
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lA<:(m)(Ai + A2 + A3) dm with 
Ai = c'{m)Qrn[Z\To > am]Qm{To > am), 

ram 

A2 = c'(m) / Q„[Z|ro = r, To = ri]Q„(ro G [r, r + dr],To = Ti), 
Jo 

ram 

A3 = c'(m) / Q^[Z|ro = r, To = r2]Q^(ro G [r, r + dr], Tq = Ta). 

JO 



We shall assume m ^ A. 



4.2.2. Computation 0/ Ai. We have, using formula ([UD 

Ai = c'(m)Q„(ro > a„)Q^[e-/o'™--x(di)n, j^.^ > -^j 
= c'(m)e-/o'""<^'W™-*))'^* 



exp 



-2/3 / {wm{t) - c{m - t)) dt - / e-^'^^"-*) ^^Evc(d^)[l - e'^^^'-W-^^™-*))" 



r /"Cm 

= c'(m) e-/o™ <f>'i<n.-mt exp |- (it[V'Le(«^™(0) " V'Evelclrri - t))] 
= c'(m) e^™-'^™ ^'m)dt ^^,,(^™{t)) 
= c'(m-a„)e-/o'^'"'^*'^Eve(-".W), 
where we used (pO|) for the last equality to get 

(24) eC--am ^'(^W)'^* = e-['°^i-'mZ--am = ^V^l^ _ 

c'(m) 

4.2.3. Computation of IS.2. Using Section [4.1.21 we get 

exp(- e-^"(™-*) £^Eve(d^)[l - e"^"']) 
Jo 

ram 

exp(- / dt e-^'^^™-*) £7r((i£)[l - e"^"*]) 



— ^'c(m— r) /)/ —i'n* 
/(0,oo) <Po(c("i-'")) 
exp(- / dt[V'Eve(^m(i)) " i^Evei^^i^ " *))]) 



exp(- / dt[^'«(t)) - ^'(c(m - t))]) 

^ r 

(pQ^c^m — r)) 
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We deduce from Section [4.1.11 

ram 

A2 = c'(m) / Qn,[Z\TQ = r,To = Ti]Q,n{To G [r, r + dr],To = Ti), 
Jo 

= c'{m) / dr UwlXr)) e" ^0 exp(- / dtW^^,{w^{t)) - ^'^^Mm - t))]) 

Jo Jo 

ram 

exp(- / dt[i^'{w*^{t)) - V'(c(m - t))]) 

= c'{m) ^'(cim-t)) dt / ^'^(wUr)) e- ^0 '^EveC-^W)"/;™ ^'(v^Lit)) 

Jo 

ram 

= c'{m - dm) / dr ^'o{w*m{r)) e" ^* ^LeC-^^W)-/;'" 

JO 

where we used (pH) for the last equality. 

4.2.4. Computation 0/ A3. Using Section [4.1.21 we get 

Qm[Z\To = r,To= T2] = e-2/^i'o("'™(*)-c('"-t)) dt-2(3 j;"' {w^(t)-c{m-t)) dt 

exp |- (ite-^'=("-*) £7rEve(d^)[l - e"^"*] 

r ram "] 

exp | - y dte-^"^""-'^ £TT{di)[l - e"^"*] 

= exp(- / dt[V'Eve(^^m(t)) " V'Evelcl"! " *))]) 
JO 

exp(- / dt[^'{w*m{t)) - ^'{c{m - t))]). 



We deduce from Section 14.1.11 

ram 

A3 = c'(m) / Q„[Z|ro = r, To = TsIQ^ITq G [r, r + dr], Tq = Ta), 

JO 

= c'(m) / dr ai^^ e" /o ^'C'^C— *)) exp(- / dtW^,,[wm{t)) - i^'^Uc{m - t))]) 



exp(- / dtW{w*m{t)) - i^'{c{m - t))]) 

J r 

c'(m) e^"" i''i<m-t)) dt / ^^^^ g- /J- dt ^Lc("'-(*))-r'" 'i* '^'("'^W) 
Jo 

' ^ir ai^^ e- /o ^'^^^M))-fr dt ^' ^ 

where we used ([21]) for the last equality. 
4.2.5. Computation of A2 + A3. We have 

ram 

A2 + A3 = c'(m - o™) / dr ^\w*^{r))e--fo'^t^'B..(^'^m-irdti,'{wUt)) . 



c'(m - am) 



Differentiating (|2ip w.r.t. time and m, we get for t < m 

Om{w*my{t) - dmW*m{t)i^'{w*m{t)) = 0. 
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Notice also that for m > t > Um, we have dmW*{t) = c'{m — t) and thus 

dmW*{am) = c{m- Um)- 

We get 

exp(- I Mm = ^-^^ = ^J^^^y 

Differentiating (|22p w.r.t. time and m, we get for t < m 

We deduce that 

A2 + A3 = / drdmwUt)cp'iw*^{r)) e- V-LeC-- W) 

a^u;„(r)e-i'o '^*^Evc("'-(*)) 

= 5^u;™(0) - dmWmid^) e- /o''" V-LeC-^^W) . 

Notice also that for m > t > one has dm'w{t) = c'{m — t), in particular dmw{dm) = 
c'{m — dm)- This implies that 

A2 + A3 = a^«;„(0) - c'(m - am) e- ^0"" ^LeC-'-W) . 

4.3. Conclusion. Thus, for m ^ A, we have 

A1 + A2 + A3 = a„u'„(0), 

and 

5„u;„(0) = c(e) + u;oo(0) - u;,(0) = u;(0). 
This ends the proof of the Theorem. 

Acknowledgments. The authors wish to thank an anonymous referee for his numerous 
and useful comments which improved the presentation of the paper. 
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